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Chiral theories with spontaneous symmetry breaking such as the Nambu-Jona-Lasinio (NJL) 
model lead to the existence of a scalar mode. We present in a detailed manner how the corresponding 
low momentum effective lagrangian involving the scalar field can be constructed starting from the 
NJL model. We discuss the relevance of the scalar mode for the problem of the nuclear binding 
and saturation. We show that it depends on the nucleon mass origin with two extreme cases. 
If this origin is entirely due to confinement the coupling of this mode to the nucleons vanishes, 
making it irrelevant for the nuclear binding problem. If instead it is entirely due to spontaneous 
symmetry breaking it couples to the nucleons but nuclear matter collapses. It is only in the case of 
a mixed origin with spontaneous breaking that nuclear matter can be stable and reach saturation. 
We describe models of nucleon structure where this balance is achieved. We also show how chiral 
constraints and confinement modify the QCD sum rules for the mass evolution in nuclear matter. 
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I. INTRODUCTION 

The relation between the fundamental properties of low energy QCD, namely chiral symmetry and confinement, and 
the nuclear many-body problem is one of the most challenging aspect of present day nuclear physics. One question is 
how the interplay between chiral symmetry and confinement in the nucleon structure manifests itself in the nuclear 
many-body problem. In a set of recent papers [IHH we have associated the mean-field nuclear attraction with the 
in-medium modification of a (chiral invariant) background scalar field which reflects part of the evolution of the chiral 
quark condensate. In this framework the nuclear medium can be seen as a shifted QCD vacuum. Nuclear stability is 
ensured with the phenomenological incorporation of the nucleon response to this scalar field. This response depends 
on the quark confinement mechanism inside the nucleon. This framework has been implemented in nuclear matter 
calculation at the Hartree level In a subsequent work [2[ we also incorporated non relativistically the pion loop 
correlation energy. A full relativistic Hartree-Fock (RHF) calculation was then done in Q allowing to reproduce also 
the asymmetry properties of nuclear matter. 

The aim of this paper is to discuss the foundations of this picture and the nature of this scalar background field. 
Although the concept of a scalar field has been widely used for nuclear matter studies [|| its precise origin or meaning 
is still a controversial subject. The problem is that there is no sharp scalar resonance which would lead to a simple 
scalar particle exchange. In our approach instead we stress the chiral aspect of the problem. As soon as we start 
from a model which gives a correct description of chiral symmetry breaking in the QCD vacuum such as the Nambu- 
Jona-Lasinio model (NJL), the emergence of a scalar field linked to the quark condensate cannot be avoided. This 
is by construction a low momentum concept which does not imply the existence of a sharp scalar meson if the effect 
of confinement is taken into account. Indeed it has been demonstrated by Celenza et al [a, H| that the inclusion of a 
confining interaction on top of the NJL model pushes the qq scalar state, located originally at twice the constituent 
quark mass, well above one GeV. As for the /o(600) it appears as a broad resonance 7T7t resonance which has no direct 
relation with the scalar field. The explicit construction of the scalar field can be done using a bozonization technique 
based on a derivative expansion valid et low (space-like) momenta. The corresponding "scalar mass" , which is around 
twice the constituent quark mass, is in a fact a low momentum parameter related to the inverse of the vacuum scalar 
susceptibility. We remind for completeness that, according to [|1 0, the confining interaction has little influence on 
the low momentum parameters entering the effective lagrangian. 

A priori the range of this mass and the magnitude of the scalar coupling to the nucleon make it relevant for nuclear 
physics. The real question for this relevance is intimately related to the problem of the structure of the nucleon and 
the origin of its mass. The respective roles of spontaneous symmetry breaking and confinement in the generation 
of this mass are indeed crucial. Confinement has little effect on the low momentum parameters but it leads to the 
concept of a nucleonic response to the scalar field, as was originally introduced by P. Guichon Q. Without it, (i.e., 
in the pure NJL model), nuclear matter would not be stable and would collapse Q, due to attractive three-body 



2 



forces (tadpole diagram). On the other hand if the nucleon mass were entirely due to confinement, as in the MIT 
bag model, the background scalar field of the NJL model would be irrelevant for nuclear physics since its coupling to 
the nucleon would vanish. The reason is very simple as the quarks inside the bag, a bubble of perturbative vacuum, 
do not feel the presence of the surrounding scalar field. Said differently the constituent quarks to which this scalar 
field couples are in this case absent in the nucleon. It is likely that the nucleon mass has a mixed origin, in part from 
chiral symmetry breaking and in part from confinement. In this case, the nucleon mass in the nuclear medium can 
feel the presence of the scalar field of the NJL model. At the same time it reacts against this field and it is possible to 
stabilize nuclear matter. Nuclear saturation may then result from a delicate balance between the influence of chiral 
symmetry breaking and confinement in nucleon structure. It is clear however that the importance of the role played 
by the background scalar field in the nuclear binding and saturation should not be left to prejudice and beliefs but to 
facts which may help elucidate this role. This is among the purposes of this article. Some questions to be answered 
are : if the scalar field is an actor in the nuclear binding and saturation problem, is this role quantitatively compatible 
with nuclear phenomenology? What informations do we have on the role played by confinement? Is it compatible 
with acceptable models of the nucleon structure? The last question is a motivation for the second part of this work 
where we propose models of the nucleon where confinement and chiral symmetry breaking contribute roughly equally 
to the mass. The influence of the spontaneous breaking of chiral symmetry is large enough for the background scalar 
field to act as a source of nuclear attraction. But the confining aspect is sufficient to stabilize nuclear matter. 

Our article is organized as follows. The second section is devoted to the Nambu-Jona-Lasinio model. After a 
brief reminder of the basic properties of the model, we derive from it an effective lagrangian which is valid for low 
(space-like) momenta relevant for nuclear physics studies. We first use a sharp momentum cutoff and in a final step 
for practical calculations we use a delocalized version. We conclude this section by some general comments concerning 
in particular the evolution of the nucleon mass and quark condensate when confinement effects at the level of the 
nucleon structure are incorporated. In the last section, based on a simple quark-diquark string model, we discuss how 
the interplay between chiral symmetry breaking and confinement in the nucleon structure influences nuclear matter 
binding properties. We also discuss the influence of the modeling of the confinement mechanism inside the nucleon. 



II. THE SCALAR BACKGROUND FIELD FROM THE NAMBU-JONA-LASINIO MODEL 

A. The standard NJL model 



We first introduce the NJL model in the light quark sector whose original aim is to describe the low mass mesons: 
the pion, the sigma, the rho, the a\ and the omega mesons. The lagrangian is : 



ip (i 7^9^ — m) ip + 
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(ipip) 2 + (ipij 5 Tip) 2 
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(ip^fip) 2 + (ip Y^bTip) 2 + (ipY ip) 2 



(1) 



Using path integral technics it can be equivalently written in a semi-bozonized form 
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Here ip represents an isodoublct of quark fields, X is a scalar-isoscalar field, the matrix P = f ■ P = TjPj describes a 
pseudo-scalar isovector field. The matrix = £1 M + TjVj* contains an isoscalar (f2) and an isovector (Vj) vector fields 
and A^ = Tj ■ A\j is an isovector axial-vector field. The current quark mass is m and G\ and G2 are two (positive) 

coupling constants. P is the chiral partner of the E, Aj is the chiral partner of Vj and in the limit of vanishing m 
(chiral limit) this lagrangian is chiral invariant. Coupling to left (£ M ) and right (7?/ 1 ) electroweak currents is included 
through the replacement : 



V„ + 75 A^ -> V„ + 75^ + V M + 7s-4 M with 
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In the mean field approximation the constituent quark mass Mo is solution of the gap equation 

I f A d 3 k 
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where N c = 3 and Nf = 2 are the number of colors and flavors and Ek = y/k 2 + Mj 2 . The quark condensate (per 
flavor) is (qq) — — 4 N c Mo I\ . The second form of the I\ integral has been obtained through the introduction of 
a sharp three- momentum cutoff A. This sharp non covariant cutoff is first taken for simplicity. We will use later 
a delocalized version of the NJL which corresponds to a softer cutoff procedure. The mesons can be generated as 
collective qq modes either by applying standard RPA to the original lagrangian or by performing a second order 
expansion in the fluctuating fields of the bosonized effective action. We list here some results and for that purpose 
we introduce the integral : 



I(lu) = 2N c N f I 2 {uj) = 2N c N f 



d 3 k 



1 



(2*)* E k (4M - w*y 
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The qq scattering amplitude at zero CM momentum, A4^'(u)), in the pion channel is obtained from the polarization 
bubble Qp S (w) in the pseudoscalar channel incorporating the ir — a\ mixing: 
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The difference between I(lu) = I(u>) / [f + 4 Mq Gi I(oj)] and I(oj) comes from the mixing effect. From this result we 
deduce the physical pion mass, the pion-quark coupling constant and the pion decay constant : 



Here the factor rv = 



m l= M %( m n) 9l qq 



ml ( gi \ 
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ft = Ml I(m„) 
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which in practice is very close to one, is the residue at the pion pole. 



B. Effective theory for low- momentum nuclear physics 

The meson spectrum (scalar and vector mesons) can in principle be obtained in the previous scheme. This approach 
is notoriously unsatisfactory due to the lack of confinement: in particular unphysical decay channels of vector mesons 
in qq pairs may appear but, as discussed in the introduction we aim to derive an effective low momentum theory ( i.e. 
for low space-like momenta relevant in nuclear physics) and not to discuss the on-shell properties of scalar and vector 
mesons, in particular their physical masses. Hence our resulting mass parameters for scalar and vector mesons will 
not be the on-shell masses but simply mass parameters associated with the inverse of the corresponding correlators 
taken at zero momentum. As emphasized by Celenza et al 0, 0] confinement, which is needed to prevent unphysical 
decays of mesons, plays a minor role for the low momentum fields (in particular the scalar one) relevant in nuclear 
physics. We now describe the technical steps needed to establish the form of the effective low-momentum lagrangian. 

a. Effective lagrangian from NJL model. The aim is to establish a low momentum lagrangian in the meson sector 
so as to generate the dynamics of the scalar field. Technically this can be done by integrating out quarks in the Dirac 
sea using a path integral formalism. The physical meaning is simply a projection of qq vacuum fluctuations onto 
mesonic degrees of freedom. 

In the spirit of our previous works [J-H, Q we first go from the cartesian representation, (£, P), to a polar repre- 
sentation, (<S,7r), by making the change of variables : 

m + Z + iP = SU with U = ( 2 =e i ^. (7) 

It is convenient to introduce a new quark field q defined by : 

g = £ 5 V>, q = ^k with £ 5 = e *™W2 (8) 

and new vector and axial vector fields according to : 



V* 1 



(9) 
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with 

vi 1 = \ (we + e 1 ^) , a% = \ (we - em) . (io) 

With these new fields the semi-bozonized lagrangian takes the form : 

C = q[i-fdp - S - r (V„ + 75 AJ] q-J^-trj (S 2 - mS (U + [/+)) 

+ -^r trf (fy» + V? - V^) 2 + (A 1 * + A» - A^) 2 ) . (11) 

The next step is to integrate out the quarks in the Dirac sea. In that way the kinetic energy term of the mesons fields 
will be dynamically generated from the quarks loops, i.e, from quantum fluctuations. For convenience we go from 
Minkovski space to Euclidean space. Using standard transformation rules the corresponding Euclidean lagrangian is : 

L E = qDq with D = irf ■ II M + S, 11" = P M - r„ = -id? - (Vn + 75^)- (i 2 ) 
The Euclidean partition function is expressed in term of the fermion determinant according to : 

Z = e- SF = ( dqdqe- dixiiDq . (13) 



Ignoring its imaginary part, the effective action can be written as : 

S F = —TrlnD = — Trln(DD^) 

= ~Trln(ll 2 + ia tlv W tlv + S 2 - j^V^S) (14) 

with 

II M = P M - r M = -id^ - (V„ + l5 A^) , V^S = d„S + 2i lb A^ 

= d^ v - d v T^ - i [r M , r„] (15) 

We perform a derivative expansion valid at low momentum of the fermion determinant to second order in the deriva- 
tives. The difficulty lies in the fact that we do not make an expansion around a constant (vacuum expectation value 
of the scalar field) but we want to have a formal expansion with the scalar objet S keeping its field status, so as to 
include itspossible modification in the nuclear environment. For that purpose we use the elegant method proposed 
by Chan [lOj. The starting point is the following trick which uses translational invariance in momentum space: 

Tr [A (n„, G(X))] = Tr [e ik ' x A G{X)) e"**] = TrA (IL, + G(X)) (16) 

where A represents any operator depending of the position operator X and of the generalized momentum II. Hence 
the arbitrary four-momentum, k, can be averaged : 

Tr [A (U„ G(X))} = J 0fTrA (U, + k„ G(X)) . (17) 

It follows that the quark determinant can be calculated as : 

Sf = -\mW)l (04 Tr NG- 1 ) + ^(i + Gn 2 + G(2fc.n + a )] 

with G = (fc 2 + S 2 ) 1 , a = ia^W^u - ^V^S. (18) 

As pointed out by Chan, the introduction of the momentum integration does not disturb the full trace operation 
and offers the freedom needed for manipulations, such as cyclic permutations of the operators or integrations by part 
under the condition that there is an implicit rcgularization scheme. The essential point is that the final form for the 
action is entirely expressible in terms of covariant derivatives [11^,5] and [11^,11^] = W^, as it should. Once this is 



5 



done the explicit trace over \x) states can be performed producing the 8^(0) compensating the one coming from the 
average procedure. To perform a second order derivative expansion corresponds in practice to make a fourth order 
expansion in 11^ and second order in a. The result is : 



1 fd 4 kd 4 
Sf = ~2 



f d 4 kd A x ( ,„_ lN 1 „ 2 2 

J-^r tr DFC (ln{G ')-- 2 G-<? 



<2 r<2 C 2 ~ 



V G S z d^Sd^S + — W^W^ ) (19) 



where the trace operation acts in Dirac, flavor and color spaces. The resulting effective action is reducible to a local 
lagrangian. Coming back to Minkovski space but keeping the momentum k explicitly in Euclidean space this local 
Lagrangian, C mes , has the form : 

C mes = ± 2N c N f I 2S (S)d^S d^S - W(S) 

+ T^- tr f( U + C/t - 2 ) + ^N c N f I 2 (S)4A^-A,S 2 

4 Lrl Z 

tr/((V + V£ - V?) 2 + (A" + A£ - A^f) 



4G 
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- i 2iV c JV> hv{S) (n» v % v + W v ■ Vfiv + A» v ■ A^ . (20) 
The chiral effective potential, W(S), is : 

W(S) = -2N c N f I (S) + {S 2 ™ ] ■ (21) 
the quantity, —Jo (<S), represents the vacuum energy density per degrees of freedom associated with the Dirac sea 

M d 3 k 



h (S) = J In {kl +S 2 )=J o ^ E k (S) E k ee E k {S) = V^Ts^ 



(22) 



where the second form corresponds to the sharp non covariant cutoff. The integral 1% (S) is the usual NJL loop 
integral : 

f id*k E 1 f A d*k 1_ 

l2 ^> J (2n)± {k 2 E+S 2f~J Q (27r)34£«(5)- ^ 

The integrals I2S.V (<->) entering the scalar and vector kinetic energy terms are : 

hs = h- 2S 2 h + 25 4 / 4 , hv =h + S 2 h - \s A h 



id i k E 1 f A d 3 k 3 



(2^) 4 (fc| +5 2 ) 3 Jo (2 tt) 3 16^ (<S ) 
ld 4 k E 1 f A d 3 k 5 



4( )= i (^(fclTW^o (2tt)3 96^(5)- (24) 

This expression of the lagrangian shows that our objective of eliminating the vacuum qq fluctuations in terms of 
"observable" background fields is realized. Contrary to usual methods the integrals Iis,v appearing in the scalar 
and vector kinetic energy terms differ from the usual I 2 integral. However the difference between l2s,v an d I2 has 
no influence on nuclear matter calculation, at least in the Hartree approximation, since the derivative terms play no 
role. It is also important to notice that the derivation has been done assuming implicitly a covariant regularization 
procedure and may not be strictly valid for the non covariant cutoff which will be used in practice. 

The meson-like lagrangian written above in cq. (|20|) has been obtained by integrating out the fluctuating quark 
fields (quarks in the Dirac sea). It remains to enlarge this lagrangian to the "classical" quark fields corresponding 
to the valence quark sector. Formally this can be done by introducing quark source terms and splitting the quark 
field into a classical part Q and a fluctuating part. The integration of the fluctuation part produces the above quark 
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determinant and the lagrangian C mes . After a Lcgendre transformation, we arrive at the valence quark effective 
lagrangian which simply reads : 

Cval =Q[i - S - 7 M (V? + 75 Q. (25) 

Remind that S is still a field and its vacuum expectation value (corresponding to the minimum of W(S)) is the 
vacuum constituent quark mass, Mo, solution of the gap equation In nuclear matter its expectation value, S, is 
the solution of an in-medium gap equation modified by the presence of the nucleonic scalar density. It coincides with 
the in- medium modified constituent quark mass M (see next subsection). Its fluctuation enters the scalar exchange 
Fock term according to the treatment given in Q . In the following we will replace in eq. (f20|) the scalar field by its 
expectation value in all derivatives terms and in the extra mass term, I 2 (S) ■ S 2 , for the axial vector meson. In 
order to prepare the identification of canonical modes, we rewrite (omitting electroweak fields) the mesonic lagrangian 
as 

Cmes = ^2N c N f hs(S)d»Sd^S-W{S) 



f^Mf +C* - 2) + K5)S 2 tr,(A»-A 

1 9 1 



, tr f (V + V?) 2 - - 2N c N f I 2V (S) (O^Q^ + V» v ■ % v + A» v ■ 1^) (26) 

where we have introduced the quantities : 

I(S) = 2N c N f I 2 (S) , 7(5) = 2^ 7 2 {§) = i - . 

We now redefine the axial- vector meson field in order to eliminate the tt — ai mixing. For this purpose we introduce 
the canonical axial vector field, o^, defined according to : 

A^ + A ^ ,- N = g v a" with gl = 3/2 - (27) 

1 + 4G 2 S 2 I(S) yv yv 2N c N f I 2V (S) V ' 

where gv is the quark- vector coupling constant. Similarly the canonical vector (uj^v^) and scalar (S), fields are 
defined as : 

^=jyi/, V t * = g v v» 

S = g 0S S with 9 2 os= 2NcNf 1 l2s(Mo y (28) 

Here the quark-scalar coupling constant, gos, is defined at the vacuum point. Omitting the pirn coupling terms, the 
low-momentum effective lagrangian takes the form: 

C mes = \ * 28( f} d»Sd»S + W(S = g 0S S) 

+ ji[ trf{U + C/+ ~ 2) + ~ 1 ^ §2tr i W- A %) 
+ ^M V (cj^ + fl" • + X -M\ (3" ■ n„) 

1 (w^w^ + v" v ■ + a» v -a^). (29) 



The vector and axial- vector low momentum mass parameters are given by 

G 2 ' My 



M v = & ^ = l + AG 2 I{S)S 2 . (30) 



One defines the canonical pion field, <f> = f ■ $, through U = exp(i$/F) where the constant F is given by F 2 = 
7 (Mo) Mq . Coming back to the previous form of the lagrangian (eq. [2U)| . a direct inspection of the coupling of d^Q to 
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the axial weak current allows to identify the pion decay constant F„ with the parameter F. The vacuum pion mass 
parameter is finally obtained as : 

In nuclear matter the explicit S factor in front of the pion mass term renormalizes the pion mass. As in ref. Q 
we do not consider this effect since it is almost completely compensated by other contributions and the pion mass is 
expected to remain stable in nuclear matter. As for the (canonical) vacuum scalar mass, it is : 



d 2 W \ J 2 (M ) / 2 (Mf 



^=4 ^ =7^FT ^'o+ Mi). (32) 



We stress again that the quantities My, Ma, M Tl F t are not on-shell constants but low momentum effective lagrangian 
parameters, (i.e., taken at zero momentum). In practice however M v and F v differ little from the physical pion mass 
and pion decay constant. On the contrary the vector and axial vector mass parameters have a priori no reason to 
coincide with the physical omega, rho and a\ meson masses. We now develop an alternative approach which relaxes 
the sharp cutoff procedure. 

b. Delocalized NJL model. As seen before when a vector interaction term is added the tt — a\ mixing has the 
effect of decreasing the pion decay constant and it is not easy with the sharp cutoff to reach a sufficiently large value 
for Fn . This is one motivation to adopt for practical phcnomcnological calculations another version of the NJL model 
with a smooth cutoff function. As discussed below there is also physical motivation for such a smooth regularization 
associated with non localities. 

The non local version of the NJL model is obtained, for any channel, with the replacement : 

{^jYjif) (x) — > Jj(x) = J d i x\ d 4 X2 F c (x\ — x) F c (x — x%) -^(x^T jip(x2). (33) 

The presence of the form factor automatically provides a regularization procedure. Moreover such a delocalized 
lagrangian possesses a physical basis in terms of quark- instanton interaction (Tlj . We define the Fourier transform, 
f{p), of the form factor appearing in the delocalized currents Jj(x), 

F c (x) = J ^ e -*"/(p), (34) 

with /(0) = 1 in such a way that F c satisfies the normalization condition J d 4 x F c (x) = 1. As can be checked this 
procedure maintains the chiral invariance of the interaction. For actual calculations we will take the non covariant 
version of the non local NJL model : 

F c (x) = S(t) F(r), f{p) = f(p)= [ d 3 r e~^ F(r). (35) 



One practical consequence is that the interaction when written in momentum space is modified according to : 

GjS® (pr + p 2 - P3 - Pi) G 3 S^ (fx + p 2 -PZ- Pi) f(pi)f(P2)f(p 3 )f(Pi). (36) 

As it is always implicitly done in NJL (with various cutoff prescriptions), we also apply the derealization procedure 
to the current quark mass term: 

m r ip(x)ip(x) —> m J d^xi d 4 X2 F c (xi — x) F c (x — X2) ij>(xi)'i()(x2). 

It can be checked that the QCD realization of explicit chiral symmetry breaking is not affected in the sense that the 
operator identity, [Qi [Qj, H]] = H x sb, is still realized and consequently the GOR relation also holds. In the following 
we choose a gaussian for the form factor 

f(p) = e^ (37) 

where A (possibly related to the inverse of the instanton size) is the cutoff parameter of the order of 1 GeV. One 
advantage of the non local version is the smooth momentum dependence of the constituent quark mass, in agreement 
with lattice calculation. Indeed the gap equation in vacuum writes : 

M(p) = M f 2 (p), with M = m - 2 G x (qq) 
d 3 k 2 r M(k) 



(qq) = -N c N f I ^3 f 2 (k) t^L Ek = ^WTMHk). (38) 
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The dclocalized version of the semi-bozonizcd form of eq. 2 writes : 
C(x) = -ip(x)iyi J, d lx il>(x) 

d i x 1 d i x 2 4>(x 1 )F c (x 1 -x) (E + iP 7 5 + + 75^)) {x)F c {x-x 2 )^{x 2 ) 

" 4^ trf ((S " + (x) + JG 2 tTf + MA ») {X) (39) 

In the presence of the form factor we found more convenient to perform the quark integration with the original quark 
field. Going again in Euclidean space the Dirac operator, D l , is defined according to 

rl = F(p) (^ + 75 i M )F(P) 

W l = F(P)WF(P), W = T, + i~f 5 P = SU 5 ==S£l (40) 

where F(P) is an operator diagonal in momentum space whose eigenvalues coincide with the form factor f(p). We can 
see that the effect of non local coupling is to transform all the field operators 0{X) into O (X) = F(P)0(X)F(P). 
Inclusion of electroweak coupling is done by making the replacement of the type: F(P)A^F(P) — ► F(P)A fJi F(P) +A tl . 
The euclidean effective action is 

S F = —TrlnD 1 = —-Trln(D l D 1 ^) 

= - X -Trln (rf 2 + ia^W 1 ^ + W l W 1 ^ - J^W 1 ) (41) 

with 

v\ L w l = d^ 1 - [a\ l ,p 1 ) + l7 5 (d^p 1 - [v^,p 1 ] + 

= 8^1-8^1-1^1}. (42) 

For the derivative expansion we again use the momentum averaging trick. It involves F(k + P) terms and consequently 
term with derivatives of f(p). Here we take the prescription of neglecting them. Hence within this approximation all 
the fields such as S will be multiplied by the number f 2 (k) in the various momentum integrals. The calculation of 
the quark determinant is formally identical to the previous case : 

Sf = ~J ^0tr[ln(G- 1 ) + ln(l + GU 2 + G(2k-U + a t )] 

with G= (k 2 + f 2 (k)S 2 y\ a^ia^W^-j^W 1 ^. (43) 

Coming back to Minkowski space we obtain a new effective Lagrangian. We omit again in its expression pirn terms 
and terms involving couplings of pion and axial- vector fields to the derivatives of the scalar field : 



£Ls = \2N c Nf$ s (S)di>Sd ll S - W l (S) + ^-tr f (U + - 2) 



4Gi 



+ ^2N c tr f [l l 2 (S)S 2 8^U 8^ + 4/^(5) d^PA^ + AI l 28 (S) i^S 2 ) 

+ 7TT trt<y»% + A^Afj) - \ 2N c N f I 2V {S) (& v % v + V» v ■ + A» v ■ A 

4 Lr 2 D V 



(44) 



The new chiral effective potential, W l (S), is 



W l (S) = -2N c N f I l (S) + {S 2 ™> . (45) 
The quantity, —Iq(S), represents the vacuum energy density associated with the Dirac sea : 

4 = J J^i ^ {k% + f 4 (k)S 2 ) = J ^ E k (S) E k ee E k (S) = V^ + TW- (46) 
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where the second form corresponds to a non covariant prescription. The various other integrals are : 



1 2S 



Jl 



I 1 - 
1 3 — 



J 1 - 
J 4 — 



J 1 - 
1 26 ~ 



O c2 W i 
— Zo i 3 + 


O C 4 tI tI 
ZD 1 4 , 1 2V 


jl 1 c2 jl ^ o4 tI 


d A k E 


f\k) 


/•°° d 3 fc / 4 (fc) 


( 27r ) 4 


+ f 4 (k)S 2 f ~ 


/o (2 tt) 3 4^(5) 


d A k E 


f\k) 


/•°° d 3 /c 3/ 4 (fc) 


( 27r ) 4 (k% 


+ / 4 (fc)5 2 ) 3 


Jo (2 71)3 16^(5) 


d 4 k E 


f\k) 


/■°° d 3 fc 5/ 4 (fc) 


( 27r ) 4 (k% 


+ / 4 (fc)5 2 ) 4 


7 (2 71)3 96^(5) 


d A k E 


f 6 (k) 


f 00 d 3 fc / 6 (fc) 


( 2 ^) 4 (k% + nk)S*f ~ 


J (2tt)3 4£g(5) 


d 4 k E 


f s (k) 


f 00 d 3 fc 


( 2 ^) 4 {k% 


+ / 4 (fc)5 2 ) 2 


7 (27r)34£g(5)- 



28 7 (270^1 + /4(jfe)52) 2 7o (27r)3 4£g(5)- 1 j 

Again the pion-axial-vector mixing has to be eliminated. Keeping only the relevant terms for nuclear physics purpose 
we obtain the following effective lagrangian 



-l ± ± 2S 

"mes T i 



d^SduS - W{S = g 0S S) 



1 4s (g) 

2 4s ( M o) 

+ jF 2 M^ + [/t _ 2 ) + ^/(^^Wl 

+ ^4 (uA^ + if ■ ^) + ^Ml (3" • a M ) 

- i (w""^ + «^ • <V + 3"" • a^) (48) 
where we have introduced the quantities 

I 1 (S) = 2N c N f 4 (5) , P (S) = 2N c N f P 2 (5) 

Af c2 jl2( c\ 

P(S) = 2N c N f 4(5) = j'(S) - 1+ * G J 2 l i iS y F " = P ( M o) M l (49) 

The constant F given above defines the canonical pion field $ through U = exp(Q/F). The canonical vector and 
axial vector fields are defined according to 

& = g v ^ V»=g v v» A" + * G . 2 f L = 9v a» (50) 

1 + 4 G 2 6^4(6) 



with 



3/2 2 = sf, Mj 

2N c N f P 2V (S)' v G 2 ' M 2 



g^ = ^7^ FoT ' ^ 2 =^, ^ = 1 + 4G 2 4(5)5 2 . (51) 



Coming back to the previous form of the lagrangian (eq. |4"4"]) . the axial weak current, introduced through the 
replacement / 2 (fc)A /J — > / 2 (fc)A Al + Aa, coupling to the quantity 9^$ allows the identification of the pion decay 
constant parameter F^ with the parameter F: 



P (Mq) M 2 _ 2_4(Mp)Mg 2G 2 M 2 4(M ) 

F F 1 + 4G 2 Af 2 4(M ) " 

The pion mass parameter keeps its formal expression : 
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As for the scalar field the rcscaling parameter (similar to eg . (|28l) ) becomes g 2 s = (2N c Nf /^(Mo)) 1 and the 
corresponding canonical mass parameter is : 

Notice that the results differ from the sharp cutoff case since we now have three different integrals l\ , I 26 , I 2S in place 
of one I 2 ■ 

We have a priori four parameters, G\, G 2l A and the bare quark mass m. We use 

A = lGeV, m = 3.5MeV, G\ = 7.8 GeV~ 2 . 
and we obtain for the vacuum quark mass at zero momentum and the quark condensate : 

M = 371 MeV (qq) = -(286 MeV) 3 . 

The G 2 parameter constrained to reproduce to reproduce the VDM phenomenology is : 

\Vdm _ 9y _ ( 2-65 x 
Ml V 0.770 



With this value we obtain for the pion parameters : 

G 2 = (G 2 ) VDM => = 91.3MeV, M n = 141.3 MeV. 

However in nuclear matter calculations we allow for a small variation of G 2 : 

G 2 = 0.78 {G 2 ) VDM F v = 93.6MeV, M v = 137.8 MeV. 

With this set of values the low momentum mass parameter for M a defined in eq. (|54[) and for other quantities are : 

M a = 923 MeV, M v = 1256 MeV, M A = 1398 MeV. 

The numerical value of the vacuum scalar coupling constant is gos = 5.55. Notice that the numerical values of these 
mass parameters and the associated coupling constants are significantly altered by the fact that hsy differ from I 2 . 
However the ratios gv/My and gos/M a = 0.006 MeV -1 which are the relevant quantity for nuclear matter calculation 
are not sensitive to this effect. 

It is interesting to derive the potential of the equivalent linear sigma model. It is obtained through a second order 
expansion in S 2 of the Dirac sea energy defined in eq. (|4"51 14"6")) around its vacuum expectation value M5 2 . We recover 
the usual linear sigma potential once we introduce a rescaled '"effective"' scalar field (S) e ff = {F 7T /M )S, normalized 
to Ftt in the vacuum : 



W L 



crM 1 {MDeff ~ Ml \ _ (M 2 ) eff -Ull 

(b Uff K (M2) e// - Ml 



4 2F1 



F n M*(S) eff . (55) 



This potential has the form of the linear sigma model potential and the parameters, instead of being arbitrary 
constants, have been dynamically generated. We remind that our field (S) e ff is a chiral invariant, so as to respect 
all chiral constraints for the mass evolution, while the scalar field of the sigma model, a is not. The mass associated 
with this effective scalar field is : 

(M 2 ) eff = U M 2 ( A j + Ml J {M a ) eff = 659 MeV 

\ \ / vac / 

and the corresponding scalar coupling constant is the one of the quark level linear sigma model (gos)eff = Mq/F v . 
Taking the nucleon as a naive juxtaposition of three constituent quarks, its mass evolution at low density goes as 
follows : 

M* N ~ M N + 3(g s)eff {(S) eff - F n ). 



11 



c. Completion of the model. In order to prepare the ground for a quark-diquark model of the nucleon (see next 
section), we also introduce an interaction in the quark-quark channel. We limit ourselves to the color 3, scalar- isoscalar 
diquark channel. The corresponding contribution to the NJL interaction is 

Q 

Cdiquark = — (lpcij5T2/3ail) (if) i^5T 2 /3 a Ipc) (56) 

where /3 a = y/3/2X a (a = 2, 5, 7) are color matrices, i\) c = 172V 1 * is the charge conjugate of the quark spinor and 
G\ is a (positive) coupling constant. It can be generalized to the delocalized version exactly as for the case of the 
interaction in the qq channel. In the presence of the diquark channel, the bozonization procedure can be also done 
in presence of diquarks using the Nambu-Gorkov formalism. In the simplest approximation with a constant scalar 
background field, one obtains the mass and kinetic energy lagrangian for the scalar-isoscalar diquark fields, Ag : 

C l diquark = 3"A§ d,A a ; - Ml A a s Af. (57) 

The diquark mass in a background scalar field S is: 

Ml (S) = \— - 2 4^T- (58) 

We now come come to the valence quark sector of the lagrangian. It also has an explicit delocalized form written 
below : 

£ l val (x) = ^{x)i-fd^{x) 

d A x x d A X2^{x x )F c {xx-x) (S + iPj 5 + + 75 (x)F c (x-x 2 )V(x 2 ) 

d 4 x l d i x 2 F c (xi - x) i lb r 2 f3 a A$(x) F c (x - x 2 ) *(s a ) 

d 4 Xl d 4 x 2 *(.ti) F c ( Xi - x) i l5 r 2 p a A a s (x) F c (x - x 2 ) * e (x 2 ). (59) 

The full lagrangian, C l mes + C l diquark + C l val , can be utilized to describe the nucleon, generating models from a simple 
juxtaposition of constituent quarks to more refined ones such as chiral solitons (nucleon bound by the chiral fields) or 
quark-diquark models including quark exchange diagrams. However with these models nuclear matter still remains 
unstable. This is the motivation for the next section where we introduce on top of the effective NJL lagrangian some 

confining interaction between quarks or between quark and diquark. If it is done the momentum of a valence quark 

1/2 

inside the nucleon will be limited topw K J ring ~ ^qcd ~ 200 MeV which is much smaller than the scale A sa 1 GeV 
entering the form factor. In such a case the delocalization effect is essentially not visible, and one can ignore the effect 
of the form factor on the valence quark dynamics. If wc again perform the chiral transformation on quark fields and 
vector fields, the valence quark sector lagrangian used in nucleon structure calculation can be safely taken as: 

C val ~ Q [i^ -S-'f (V„ + 75 A J] Q 

- Qc *l5T 2 p a A a J Q - Q c i l5 r 2 p a A a s Q. (60) 



C. Concluding remarks on this section 



In the pure NJL picture, at finite baryonic density, the value of the constituent quark mass, which is the expectation 
value of the scalar field 5, is modified. It can be obtained self-consistcntly from a gap equation modified by the 
presence of a Fermi sea. However in the real world baryonic matter is not made of independent constituent quarks 
but of clustered objects, the nucleons. These nucleons are embedded in the scalar background field, S, and the nuclear 
medium can be seen a priori as a shifted vacuum. The nucleon mass will depend in some way on the scalar background 
field and the energy density of symmetric nuclear matter at the Hartree level reads : 

V =£ " = j &iPF _ P) (\/p 2+M n& - ( M nU?J + W{S) + 9 ^ P \ (61) 
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The expectation value for the scalar field is self-consistcntly obtained by minimization of the energy density, 

|j=0 * ^ = -2(m)(S)-^(S)pf(M N (S)), (62) 

which constitutes an in-medium modified gap equation. The connection between the field 5, normalized to the quark 
mass, and the scalar field s used in our previous work [Ml is s = (Fir / M ) (S- M Q ) and p?(M N (S)) is the nucleonic 
scalar density. The scalar coupling constant of the nucleoli to the effective scalar field (which is normalized to F„ in 
the vacuum) is : 



. Mo (dM N \ 



which depends crucially on nucleon structure. For instance if the nucleon mass fully originates from confinement (bag 
models), ( d g I £ ) = 0, the scalar field just decouples from the nucleon, {{gs)eff = 0). In this case there is no shift of 
the vacuum and the scalar field is thus an irrelevant concept for nuclear matter studies. On the other extreme if the 
nucleon mass fully originates from chiral symmetry breaking (naive additive NJL, chiral soliton), then the nucleon 
mass in the medium is affected by the scalar field associated with the dropping of the chiral condensate. However 
in this case the attractive tadpole destroys stability. Only in the case where the nucleon mass has a mixed origin, 
the scalar background field can contribute to the nuclear attraction without destroying the stability and saturation 
properties. In that case by rearranging its quark structure linked to the confinement mechanism, the nucleon reacts 
against the scalar field generating effectively repulsive three-body forces. The origin of this repulsion lies in the decrease 
of the scalar coupling constant of the nucleon. In short a possibly important part of the saturation mechanism is 
associated with the progressive decoupling of the nucleon from the scalar field associated with the dropping of the 
chiral condensate. In the next section we will introduce nucleon models capable of achieving the balance between 
large enough attraction and sufficient reaction. Of course, one falls here in the modeling uncertainties. However we 
show below that a stringent constraint exists for the numerical value of the scalar nucleon coupling constant which is 
model dependent, from the value of the free nucleon sigma commutator. 

The pion-nucleon sigma term is an important piece of experimental information. It is obtained from the Feynman- 
Hcllman theorem : = m (3Mn /dm) ~ 50 MeV. It receives a contribution from the pion cloud, a^ lon cloud \ 
According to previous works @, [H|-[l4| we expect: a^ lon cloud "> ~ 20 MeV which corresponds to a pion cloud self-energy 
of —420 MeV . For the non pionic part an explicit calculation in the NJL model shows that the linear sigma model 
result is recovered but with the nucleon structure aspect hidden in the scalar coupling constant, g a = (gs)eff(Mo) : 

(nopion) _ p Mg 

- F * 9 °mu~f- (63) 

Its numerical value has to be a ^:°P lon ^ — aN _ (x& lon cloud > ~ 50 — 20 ~ 30 MeV. This separation of the sigma term 
into two pieces is quantitatively supported by the lattice study of Leinweber et al [HJ on the nucleon mass evolution 
with the bare quark mass. In their work the pionic part of this evolution which has a non-analytical behavior is 
calculated explicitly and subtracted out. For the rest an expansion is made in powers of m\. The linear term in m\ 
is linked to the non pionic sigma commutator, giving a value cr &° pl,on > ~ 29 MeV 1 close to our value. With the value 
(M a ) e ff = 659MeV given previously it leads to g a ~ 7. 

Concerning the nucleon mass problem, there exist QCD sum rules which link in an approximate way the nucleon 
mass to the condensate both for a free nucleon [l5| and for a bound one [TH, H3] ■ For a dilute medium, these sum 
rules lead to the following mass evolution : 

((nopion) . (pion cloud)\ 
+a " ) p „ , , ± _ <£!> m) 

M N -(qq) vac 1 F*M* 1 FlMl ~ F„ 2F% [ > 

where the last expression is the one obtained from the condensate evolution in the NJL model. We see that the 
quark condensate modification receives two contributions, one from the scalar field and one from the pion cloud, 
reconstituting at low density the full pion nuclear sigma term. 

Our description brings important restrictions to this expression. Firstly the pion cloud contribution to the sigma 
commutator sigma commutator, a (P lon cloud ^ contributes to the condensate evolution. It does not contribute to the 
mass evolution, otherwise chiral constraints would be violated (such as the presence of a term in m T in the NN 
potential, forbidden [l8| by chiral symmetry) . In fact its influence on the mass vanishes in the chiral limit and hence 
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it is a small effect which we ignore. Only the scalar piece, a™^ ton , should then enter the mass evolution. In a pure 
chiral theory such as NJL the low density expansion of the mass evolution is then : 

This separation already introduces a model dependence in the prediction of the mass evolution. This expression would 
hold for instance for an assembly of nucleons described as superpositions of constituent quarks or by chiral soliton 
models. Secondly when we introduce confinement, i.e, when we go beyond NJL, the evolution of the nucleon mass 
with density, also depends on the nucleon structure through the value of the coupling constant of the scalar field to 
the nucleon, g a . Indeed, to leading order, we have = Mn + g a s with s = —g a p/ {M%) e ff which gives : 

M* J-, (nopion) (nopion) 

M N ~ M N F% Ml - 10 F2 M 2 1 ' 

In the linear sigma model where g a = Mpj/F n = 10 we recover the Ioffe sum rule generalized at finite density just 
corrected from pionic effects. With confinement the value of the scalar coupling constant is reduced and the mass 
evolution is slower than the condensate one. The suppression of the pionic contribution to the mass evolution further 
accentuates the difference between the mass and condensate evolutions. For instance at normal nuclear density the 
condensate has dropped by ~ 30%. With the value g a ~ 7 deduced above, the mass reduction is significantly lower, 
~ 13%. 



III. EFFECT OF CONFINEMENT: SIMPLE MODELS FOR THE IN-MEDIUM NUCLEON 

We now come to the last point of this paper, namely the modeling of the nucleon mass origin and the scalar response 
of the nucleon defined from the second derivative of the nucleon mass with respect to the scalar field : 



*"M-^-jl^)-¥A-^iry (67) 

For a nucleon made of the simple adjunction of three NJL constituent quarks (or a NJL quark and a NJL diquark) 
the scalar coupling constant is independent of the scalar field and there is no scalar response. The importance of 
the response is related to the respective roles of chiral symmetry breaking and confinement in the generation of the 
nucleon mass. In the following we will consider nucleons built from NJL quarks of mass M bound by some confining 
force. The information that we need is contained in the relation between the nucleon mass Mn and the NJL mass 
M, Myv = f(M). The scalar coupling constant of the nucleon is related to that of the NJL quarks, which is M/F n , 
through gs = ~§m^§§-- The next derivative with respect to s gives the nucleon scalar response. A non vanishing value 

requires 7^ and it entirely arises from confinement. For instance in soliton models where the quarks are bound 
only by chiral forces the nucleon mass is proportional to the NJL mass and the second derivative vanishes. In the 
following we establish this relation for different models showing their effect on the saturation properties. 

In a previous work [l9j we have introduced a model of a nucleon made of three constituent quarks bound together 
by a confining harmonic force. The magnitude of the scalar response which followed was too small to prevent the 
collapse of nuclear matter. We will come back later to this type of model. A possibility of improvement is to reduce 
the relative role of chiral symmetry breaking. This can be achieved by considering a nucleon made of a quark and a 
sufficiently light diquark to leave enough room for confinement. A practical advantage is that a three-body problem is 
transformed into a simpler two-body problem. Beside this simplification, there are theoretical and phenomenological 
reasons to favor a quark-diquark model of the nucleon with relatively light scalar-isoscalar diquark. For instance the 
work of Shuryak et al on hadronic current-correlation functions based on a random instanton vacuum (20| finds a 
strong attraction in the scalar-isoscalar channel leading to a diquark with a mass about 400 MeV. An axial- vector 
diquark is also found but with a much larger mass of the order of 900 MeV. It is also possible to nicely reproduce 
the light baryon spectrum [2l| while a calculation without diquark correlations predicts an abundance of missing 
resonances |22|. 

As discussed in a set of works of Bentz et al (see ref. for application to nuclear matter), it is possible from the 
NJL model to construct a nucleon with a diquark component. Introducing the standard interaction in the diquark 
channel as discussed previously one obtains for the mass of the scalar diquark result quoted in eq. (|58[) . This mass is 
also medium dependent since it depends on the constituent quark mass. Its vacuum value is strongly sensitive to the 
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FIG. 1: Mass of the quark (dashed line), of the diquark (dot-dashed line) and of the nucleon (full line) versus the relative 
deviation, $ = (Mo — M)/Mo = \s\ /-FV, of the scalar field with respect to its vacuum value. 



value of G\ which defines the quark-quark interaction (see eq. I56[) . For G\ = G\ it is exactly equal to the pion mass. 
Here we choose 

Gi = 0.92 Gi M D = 398.5 MeV 

which turns out to be nearly equal to the constituent quark mass in agreement with the work of ref. [20j . In 0], it 
was realized that to obtain a scalar susceptibilityand consequently nuclear matter saturation requires a confinement 
mechanism. An infrared cutoff [1r ~ 200MeV was thus introduced in the Schwinger proper time regularization 
scheme. Such a prescription implies that quarks cannot propagate at relative distance larger than 1/fiR, hence 
mimicking a confinement mechanism. Here we propose to incorporate confinement in a more direct way. Since the 
diquark is in an anti-triplet color state, it is physically plausible that a string develops between the quark and the 
diquark as in a QQ meson. We thus introduce a confining potential between the quark and the diquark : 

V(r) = l^Kr 2 . 

In the non relativistic limit, the problem reduces to solving the Schroedinger equation for a particle with reduced 
mass /i, placed in an harmonic potential. In this limit the mass of the (in-medium) nucleon is given by : 

M N (S) = M(S) + M D (S) + with n = 




M + M D 



We take for the string tension a standard value K = (290 MeV) 3 . We obtain for the vacuum nucleon mass 
AIn = 1304 MeV . The nucleon mass origin splits roughly into a chiral symmetry breaking component (60%) and a 
confinement component (40%). The vacuum value scalar coupling constant of this nucleon to the effective scalar field 
is g a = (gs) e ff{M ) = 7.14. This leads to the value of the non pionic piece of the sigma term : a ^° pion > = 39 MeV, as 
was required. In order to show that such a model is capable of describing the saturation properties of nuclear matter 
we calculate the energy of symmetric nuclear matter in the Hartree approximation, using eq. (|61l [6"2"]l . The resulting 
curve displays a saturation mechanism driven by the scalar nucleon response (kns, proportional to the second deriva- 
tive of the nucleon mass with respect to 5, eq. |67| which has a positive value. Said differently the scalar coupling 
constant, dMjy/dS, is a decreasing function of |s| or the density. This translates into the fact that the nucleon mass 
stabilizes or even increases with increasing |s| (sec fig. 1). However the binding is nevertheless not sufficient unless 
we decrease artificially the vector coupling constant G2 at a value much smaller than the VDM result. In order to 
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FIG. 2: Binding energy of nuclear matter versus nuclear matter density in units of normal density. The full line corresponds 
to the full result and the dashed line represents the Hartree result. The dot-dashed line corresponds to the contribution of the 
Fock term and the dotted line represents the correlation energy. All the numerical inputs are given in the text. 



improve the description, although this is not necessarily consistent with our present nucleon model, we add on top of 
the Hartree mean held result the pion loop (Fock term and correlation energy) contribution obtained in our previous 
work Q- Taking the value of G2 at the value quoted previously, G2 = 0.78 (G-2) VDM , we obtain a decent saturation 
curve shown in fig. 2. Likely a fully consistent calculation within the model of the pion loop energy would modify the 
result but a fine tuning on Gi would be presumably sufficient to recover the correct saturation curve. The lesson of 
this simple model calculation seems to confirm our previous conclusions. The confinement effect (scalar response of 
the nucleon) is able to stabilize nuclear matter and the pion loop correlation energy helps to get the correct binding 
energy. 

We have shown that an acceptable quark-diquark model of the nucleon makes plausible the role of the background 
scalar field in the nuclear binding. It is interesting to investigate if other confining mechanisms can achieve the 
same result. For this we have also studied models where the nucleon is made of three constituent quarks moving in a 
mean-field linear confining potential but shifted with a constant attractive potential mimicking short range attraction : 

V=^±^{K 2 r- 2V ). 

This model has been successfully utilized for baryon spectroscopy studies by Jena et al [23|. We do not aim to 
justify this particular equally mixed scalar and vector confining potentials, the main motivation being the existence 
of analytical solutions. The energy of the lowest orbit, solution of the Dirac equation, is : 

E(M) = M —2Vq + \/K^x q witha; g solution of x\ + 2 ~ x\ - (2.33811) 3 = 

v K2 

and the mass of the in-medium nucleon (in absence of CM correction) is Mn{§ — M) = 3E(M). Hence the quark 
mass contribution (essentially the chiral symmetry breaking contribution) to the quark orbital energy and then to the 
nucleon is reduced due to the presence of the attractive shift, —2Vq, leaving more room for the confining part. The 
scalar coupling constant (still omitting CM correction) can be written as 

( x ,^ M (dM N \ M 

where q s = J d 3 r (u 2 — u 2 ) (r) is the quark scalar charge. We see that the scalar field contribution to the sigma 
term is represented by the usual integrated scalar quark density as in bag models. In practice we also include in the 
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FIG. 3: effective scalar coupling constant versus the relative deviation, (Mo — M)/Mq = \s\ /F v , of the scalar field with respect 
to its vacuum value for the linear confining potential (dashed line), the quadratic linear potential (dotted line) and for the 
quark-diquark model (full line). 

numerical calculation the effect of CM correction using the results quoted in ref. [13] . If we take K<x = (300MeF) 2 
and Vo = 200 MeV it is possible to obtain a saturation curve but the saturation has the tendency to come too 
early. Certainly this point deserves a more detailed study. Here we wish to concentrate on the main result, namely 
a decreasing scalar coupling constant when increasing \s\ as demonstrated by the dashed curve on fig. 3. We also 
checked that replacing the linear potential by a quadratic potential, 

V = ^~y^ Q K 3 r 2 - 2 V ^j with K 3 = (300 MeV) 3 , V Q = 200 MeV, 

one obtains similar results as depicted on fig. 3 (dotted curve) . It is worthwhile to notice that this model differs from 
the one used in [l9[ by the introduction of the constant attractive shift — 2Vo- The energy of the lowest orbit is the 
solution of the equation : 

E = M -2V + -J-^^. 

2 V E + M 

Again this shift allows to reinforce the role of confinement in the origin of the nucleon mass. 

Also shown on fig. 3 is the behavior of the scalar coupling constant for the quark-diquark model. In this case, the 
decrease at low density is less strong which translates into a softer equation of state. According to a preliminary study 
based on a variational relativistic calculation the strong dropping beyond \s\ / « 0.2 (which roughly corresponds 
to normal density) might be to some extent an artifact of the non relativistic approximation. 

IV. CONCLUSION. 

We have studied the role played by the spontaneous breaking of chiral symmetry in the problem of the nuclear 
binding. The existence of a scalar field linked to the quark condensate emerges in chiral theories such as the NJL 
one. This field may be at the origin of the masses. This the case in the NJL model or the linear sigma one. In this 
case several things follow naturally. The partial restoration of chiral symmetry in dense matter implies a reduction in 
magnitude of the condensate and hence of the nuclconic mass, which could a priori account for the nuclear binding but 
in this case a tadpole term inherent in these theories destroys stability. A combination with the confining aspects is 
able to restore stability. Confinement indeed reduces the coupling constant of the scalar field to the nucleon and makes 
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it field dependent. Equivalently it introduces a scalar response a the nucleon to this held in such a way that the nuclear 
medium reacts against a build up of the scalar field with increasing density, which helps in the saturation problem. 
However confinement should not be the only origin of the nucleon mass since in this case the scalar background field 
decouples from the nucleons. It is only in a mixed case, with a simultaneous influence of spontaneous symmetry 
breaking and confinement that the scalar field can be an efficient actor in the nuclear saturation problems. We 
have given examples of nucleonic models where this balance is achieved. They require the role of confinement in the 
generation of the mass to be sufficient. We have shown how confinement affects the QCD sum rule for the in medium 
nucleon mass originally shown to follow the condensate evolution. Our formula shows that the mass evolution is 
reduced as compared to the condensate one by a factor r, ratio of the scalar coupling constants in the presence and 
in the absence of confinement. In addition, as we pointed out in previous works, only the non pionic part enters the 
mass evolution. These combined effects considerably reduce the mass evolution as compared to the condensate one. 
Nevertheless the remaining effect can be sufficient to make the scalar field of chiral symmetry breaking an important 
actor in the nuclear binding question. 
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